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INTENSITY  OPTICAL  CORRELATION  THROUGH  OPTICAL  DIFFERENTIATION 
PREPROCESSING 

CHEN  WENSHANG,  WANG  ZHAOQI,  MU  GUOGUANG 

Institute  of  Modern  Optics,  Nankai  University,  Tianjin 


A  high  discrimination  intensity  optical 
correlation  system  is  presented  in  this 
article.  By  using  a  lensless  spacial 
differentiation  filter  to  achieve  edge 
enhancement  for  the  input  images,  the 
discrimination  of  the  intensity  optical 
correlator  with  the  Fresnel  holographic  filter 
is  improved.  Theoretical  analysis  and 
experimental  demonstration  are  given. 


The  use  of  a  Fresnel  holographic  filter  (FHF)  to  achieve 
intensity  optical  correlation  discrimination  has  been  studied  in 
some  depth[  1,2f  .  This  new  type  of  wave  filter  has  a  number  of 
advantages . 

However,  using  a  FHF  intensity  optical  correlation  system 
calculates  the  correlation  of  light  intensity  factors  and  not  of 
complex  amplitude  factors.  Therefore,  it  reduces  the 
discrimination  of  the  correlation  test.  Also,  for  many  real 
objects  such  as  aircraft  and  tanks,  their  Fresnel  spectrums  all 
contain  a  fairly  large  amount  of  frequency  nulls,  and  their 
autocorrelation  distribution  is  a  fairly  large  diffused  spot.  This 
makes  it  even  more  necessary  to  solve  the  problem  of  characteristic 
recognition  differentiation  (discrimination) .  The  use  of  optical 
differential  preprocessing  to  enhance  the  edges  of  images,  because 
it  eliminates  the  frequency  nulls  in  the  center  of  the  image  and 
some  low  frequencies,  it  makes  the  higher  frequency  signals  stand 
out  and  sharpens  the  autocorrelation  peaks,  which  improves  the 
image  discrimination.  This  article  simply  combines  an  ordinary  FHF 
intensity  optical  correlation  system  and  an  optical  differentiation 
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preprocessing  system,  and  obtained  an  optical  correlation  testing 
system  which  has  both  flexibility  and  high  discrimination. 

II:  THEORETICAL  ANALYSIS 

Considering  only  the  one  dimensional  situation,  expanding  the 
real  factor  T(<»)  which  expresses  the  intensity  distribution  of  a 
characteristic  signal  to  where  it  has  cosine  factors  of  frequency 
fn  superimposed.  . 

T(w)  -  Jjj_ I„( fOrec* (^)  oos(2  ovfjs)  (1) 

In  this  equation,  Inf(n)are  the  weighting  of  the  nth  level  cosine 
component.  2L  is  the  sampling  interval.  For  the  purpose  of 
simplifying  matters,  here  we  only  consider  the  effects  of  frequency 
null  components  on  correlation  peaks.  The  above  equation  can  be 
resolved  as  follows: 

T(*)-Icr8ot(^)+S  !.(/.)  reot  (^)coa(2ar/^)  (2) 

In  this  equation,  IQ  is  the  frequency  null  component.  The 
characteristic  signal  intensity  autocorrelation  distribution  is: 

I(u)~J  T(x)T’(x+u)dx 

- jIl1' K(-&)+ 5 |t(M  (it) 

/ON 

•cos(2  «/,^)oos[2 

In  this  formula,  the  triangular  factor  A  (^X^n  the  first  item  is  one 
background  provided  by  frequency  null  components.  Compared  to 
summation  items,  it  is  a  slowly  changing  factor.  The  conditions  to 
establish  equation  (3)  are  to  ignore  the  effects  of  cross 
correlation  between  the  frequencies  in  different  space. 

We  qualitatively  analyzed  the  effects  of i^t)  on  the  sharpness 
of  the  correlation  peaks.  If  we  set  the  height  of  the  correlation 
peaks  of  the  summation  items  in  equation  (3)  at  2h,  half  intensity 
width  is  a  distance  of  h  from  the  peak,  and  forms  a  slowly 

changing  background  with  an  altitude  of  2H.  Thus,  with  a 
background,  the  distance  from  the  half  intensity  to  the  peak  value 
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point  is  h+H.  In  a  condition  corresponding  to  no  background  noise, 
the  half  intensity  point  shifts  the  H  distance  down  to  bottom  of 
the  outline,  so  the  half  intensity  width  is  increased,  and  the 
corresponding  peak  becomes  wider.  Furthermore,  the  larger  the 
frequency  null  component  (that  is,  the  higher  H  is),  the  wider  the 
corresponding  peak.  Some  low  frequency  components  in  the  image 
will  also  cause  the  same  thing.  In  order  to  obtain  a  sharper 
corresponding  peak,  it  is  necessary  to  eliminate  or  reduce  the 
frequency  null  components  and  the  low  frequency  components  in  the 
image,  and  boundary  intensity  optical  differential  preprocessing  is 
one  feasible  effective  way  of  doing  this. 

Ill:  INTENSITY  OPTICAL  CORRELATION  OF  DIFFERENTIAL  IMAGES 

Fig.  1.  Optical  Path  of  Lensless  Spacial  Differentiation  Filter 

/ 


In  order  to  achieve  intensity  optical  correlation 
discrimination  of  differential  images,  it  is  necessary  to  first  of 
all  make  a  spacial  differentiation  filter  to  use  for  optical 
differentiation  preprocessing  of  characteristic  objects  or  input 
objects  to  obtain  an  edge  enhanced  differential  image.  We  used  a 
lensless  space  differentiation  filter  as  reported  by  Qin  Qiuxiang[5] 
to  perform  the  differentiation  preprocessing.  The  lensless  space 
differentiation  filter  recording  equipment  is  shown  in  Figure  1. 
It  is  a  differentiation  filter  which  uses  double  exposure  by  one 
parallel  wave  and  one  spherical  wave  for  recording.  The  location 
of  the  focal  point  of  the  focused  spherical  wave  is  indicated  by 
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(d)  and  (a)  in  Figure  1.  The  parallel  wave  perpendicular  to  the 
hologram  H  is  exposed  first.  Prior  to  the  second  exposure,  the 
direction  of  the  parallel  wave  is  altered  by  a  minute  angle  60,  and 
the  two  exposures  occur  at  the  same  time.  With  appropriate 
developing  and  fixing,  we  obtain  the  needed  lensless  spacial 
differentiation  wave  filter.  A  two  dimensional  lens  spacial 
differentiation  wave  filter  can  be  made  using  similar  methods. 

Fig.  2:  Recording  Arrangement  of  the  Fresnel  Holographic  Filter  for 
Differential  Images 


The  differential  image  Fresnel  Holographic  filter  was  made 
using  the  optical  circuits  shown  in  Figure  2 •  The  characteristic 
object  g(x,y )  entered  through  plane  P3,  and  a  lensless  spacial 
differentiation  filter  is  placed  at  the  focal  plane  P2  behind  lens 
1^.  By  appropriately  adjusting  the  location  of  the  differentiation 
filter  (a  distance  b=A/260  along  x3,  it  is  possible  to  obtain  on 
plane  P3  an  edge  enhanced  differential  image. 

9t(x<i  yt)=*g(x,  y)*[d(x+a  +  b+>J)f,  y)-d(x+a+b,  y)] 

In  this  equation,  *  indicates  convolutional  computation,  (x,  y)  and 
(x„,  yj  are  the  coordinates  of  the  entrance  plane  P1  and  the 
differential  output  plane  P3  respectively.  Using  this  differential 
image  as  the  object  light,  and  placing  hologram  plates  on  the  P4 
plane  for  dx  behind  it,  using  the  focus  of  the  reference  light  to 
record  the  differential  image's  Fresnel  holographic  filter.  The 
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focal  point  of  the  focussed  reference  light  is  at  a  distance  of  dQ 
from  the  holographic  plate  and  a  distance  of  hQ  from  the  system's 
light  axis.  Under  Fresnel  approximation,  the  holographic  plate 
object  light  wave  and  reference  light  wave  are  obtained  as  follow: 


iit 


0(a,  £)-0i[Jsr<(®<>  <*!/< 


2J(o,  0 )-(?*«*• 


(5) 

I 

(6) 


Herein,  Cx  and  C2  are  complex  constants,  and  (a,  B)  are  the 
coordinates  of  plate  P4.  Through  appropriate  exposure  and 
developing,  it  is  possible  to  obtain  a  differentiation  filter 
Fresnel  holographic  filter.  The  amplitude  transmissivity  is 


!T(o,  0)cc  |22(o,  /S)+O(o,  /8)  |»- |S|*+ 10|9+22*0+220*  ..(7) 


In  these  four  diffractions,  one  item  which  is  related  to  optical 
correlation  computation  is 


T(a,  0) 

0)0' (a,  0) 

-ojfc*-'*** 

xjj gXx(,  d*4tyt  (8) 

—  M 


Herein,  C3  is  also  a  complex  constant. 


Fig.  3:  Intensity  Optical  Correlation  System  for  Differential 

Images 
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The  optical  circuitry  of  differential  image  intensity  optical 
correlation  systems  is  shown  in  Figure  3.  This  system  is  composed 
of  an  optical  differentiation  preprocessor  system  and  an  intensity 
optical  correlation  system.  The  object  to  be  discerned  (f(x,y) 
enters  through  plane  Px,  and  on  frequency  spectrum  plane  P2  a 
lensless  spacial  differentiation  filter  is  placed.  Through  this 
optical  differentiation  preprocessing  subsystem,  at  the  plane 
it  is  possible  to  obtain  edge  enhanced  differential  image  fM(x1_>,yK>) 

At  plane  P3  a  constantly  moving  diffusion  scattering  screen  is 
placed  to  serve  as  the  input  plane  for  the  intensity  optical 
correlation  subsystem.  At  the  same  time,  the  spacially  correlated 
light  is  converted  to  spacially  non-correlated  light.  The  Fresnel 
holographic  filter  made  as  described  above  is  placed  on  space  wave 
filter  plane  P4.  The  distance  between  plane  P4  and  P3  is  d1.  The 
intensity  optical  correlation  system  is  linear  for  the  incoming 
light  intensity,  thus  signifying  that  the  physical  quantity  of  the 
relationship  between  the  system  input  and  output  is  a  function  of 
the  irradiation  pulse  response.  We  first  calculated  the  system  s 
pulse  response  function.  Assume  that  a  single  pulse  8(<c«— y«  tf) 
enters  P3  of  the  correlation  system.  Its  light  wave  field  is  a 
single  spherical  wave  on  plane  P4.  Then  the  light  wave  field  just 
emitted  from  the  correlation  filter  plane  is 

&(«,  $)  .  (9) 

The  correlation  filter  pulse  response  factor  is  equation  (9) 
transmitted  to  the  Fresnel  diffraction  of  plane  P5: 

*(«*  !fo»  *W)-f dadfl  (10) 

— o* 

Substituting  equations  (8)  and  (9)  in  equation  (10)  and 

differentiating.  Under  conditions  d1=d0  we  can  obtain 

K®o,  yoi®',  f/)  -sr*(®o+®'— ho*  2/o+S/O  (11) 

The  irradiation  pulse  response  function  is  the  square  of  the  pulse 
response  function  module. 

hi(x o,  j/o;  I/')  -  1 0«(*o +*'“*<>>  Vo+lH  Is  (12) 
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The  system  output  irradiance  is  the  convolution  of  the  irradiance 
of  the  correlation  input  function  and  the  irradiation  pulse 
response  function:  V}  ' 

I(®o,  yo)-JJ  S/<)|  |0?(so+®<-Ao,  yo+yt)\2<faidt/i  (13) 

This,  then,  is  the  computation  of  the  intensity  correlation  of 
boundary  enhanced  differentiated  image. 

IV:  EXPERIMENTAL  DEMONSTRATION  AND  RESULTS 


Fig.  4: 


The  correlation  light  source  we  used  in  our  experiments  was  a 
60mW  He-Ne  laser.  The  recording  medium  was  a  HP  633  model 
holographic  plate  produced  by  the  Tianjin  Mircotechnology 
Development  Company.  We  used  special  "slow  developing"  processing 
technology,  and  the  diffraction  efficiency  of  this  type  of 
holographic  plate  was  30  to  40  percent.  In  our  experiment,  we 
first  used  the  picture  in  4(a)  as  the  input  object  or 
characteristic  object.  It  was  optically  differentiation 
preprocessed  through  a  lensless  spacial  differentiation  filter, 
obtaining  the  boundary  enhanced  differential  image  shown  in  4(b). 
This  differential  image  achieved  intensity  light  correlation 
discrimination  by  using  a  Fresnel  holographic  wave  filter, 
obtaining  the  autocorrelated  result  shown  in  4(c).  As  an 
illustrative  comparison,  figure  4(d)  shows  the  intensity  light 
correlation  results  of  the  original  input  image  without 
differentiation  preprocessing.  We  can  see  that  there  is  marked 
improvement  in  the  discrimination. 


In  order  to  quantitatively  explain  the  improvement  in 
discrimination,  we  also  performed  computation  of  the  simulation 
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experiment.  We  use  the  English  letters  "0",  "O’"  and  "T"  as  the 
target  object,  using  a  computer  to  simulate  the  intensity 
correlation  system  discrimination  of  the  letter  "0".  The 
discrimination  6  was  used  to  define  the  corresponding  value  of  the 
height  differential  between  the  autocorrelated  peaks  and  cross 


correlated  peaks. 


Kg 

ha ate 


Xl00% 


(14) 


In  this  equation,  Hauto  and  Hcross  refer  to  the  height  of  the 
autocorrelated  peaks  and  the  height  of  the  cross  correlated  peaks 
respectively.  In  our  experiment,  we  used  256  X  256  dot  high  speed 
Fournier  conversion  processing  of  a  128  by  128  two  dimensional 
image  (English  letters).  Without  differential  preprocessing,  the 
discrimination  between  "0"  and  "O’"  was  A=1.2%  and  for 
discrimination  between  "0"  and  "T"  it  was  0=21.0%.  However,  with 
differential  preprocessing,  the  discrimination  between  "0"  and  "0 
was  12.6  and  between  "0"  and  "T"  it  was  48.5%.  we  can  see  that 
discrimination  is  greatly  improved,  and  the  smaller  the  difference 
between  the  two  images,  as  in  "0"  and  "O'",  the  more  noticeable  the 
improvement.  Table  1  provides  the  specific  data  from  the  computer 
simulation  experiment.  The  correlation  peak  height  data  has  been 

normalized. 


Table  1:  Correlation  Discrimination  with  Original  Images  and 

Differential  images  respectively. 


Input  image 

correlation 

correlation  peak  Height 

discriminating  force 

Original  image 

Oand  C 

OfO 

1.000 

1.2 

0®C 

0.988 

Differential  image 

O' and  O' 

a%a 

0.719 

12.6 

v%c 

0.627 

Original  image 

0  and  T 

ogo 

1.000 

21.0 

0%T 

0.790 

DiffftTfttial 

0$'0'  . 

0.719 

48.5 

O' and  T* 

a%r  ; 

0.370 
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